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On the nature of the variabiUty power decay towards soft spectral 
states in X-ray binaries. Case study in Cyg X-1 

Lev Titarchuk^'^'^ and Nikolai Shaposhnikov^ 
ABSTRACT 

A characteristic feature of the Fourier Power Density Spectrum (PDS) ob- 
served from black hole X-ray binaries in low/hard and intermediate spectral 
states is a broad band-limited noise, characterized by a constant below some fre- 
quency (a "break" frequency) and a power law above this frequency. It has been 
shown that the variability of this type can be produced by the inward diffusion 
of the local driving perturbations in a bounded configuration (accretion disk or 
corona). In the framework of this model, the perturbation diffusion time to is 
related to the phenomenological break frequency, while the PDS power-law slope 
above the "break" is determined by the viscosity distribution over the configura- 
tion. The perturbation diffusion scenario explains the decay of the power of X-ray 
variability observed in a number of compact sources (containing black hole and 
neutron star) during an evolution of theses sources from low/hard to high/soft 
states. We compare the model predictions with the subset of data from Cyg X-1 
collected by the Rossi X-ray Time Explorer (RXTE). Our extensive analysis of 
the Cyg X-1 PDSs demonstrates that the observed integrated power decreases 
approximately as a square root of the characteristic frequency of the driving os- 
cillations i>dr- The RXTE observations of Cyg X-1 allow us to infer Pdr and to as 
a function of u^r- Using the inferred dependences of the integrated power of the 
driving oscillations Pdr and to on u^r we demonstrate that the power predicted 
by the model also decays as Px,diff oc i^J^'^ that is similar to the observed P^ 
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behavior. We also apply the basic parameters of observed PDSs, power-law index 
and low frequency quasiperiodic oscillations, to infer the Reynolds (Re) number 
from the observations using the method developed in our previous paper. Our 
analysis shows that Re— number increases from values about 10 in low/hard state 
to that about 70 during the high/soft state. 

Subject headings: accretion, accretion disks — black hole physics — stars:individual 
(Cyg X-1) :radiation mechanisms: nonthermal — physical data and processes 



1. Introduction 



The main goal of the presented work is to explain a decay of the emergent time vari- 
ability of X-ray emission in compact sources when these so urces evolve from low/hard (LH) to 

high/s oft (HS) states [see lRemillard fc McClintock I ( l2006l ) and lTitarchuk. Shaposhnikov fc Arefiev 
( 120071 ). hereafter TSA07, for details of observations]. In particular, power density spectrum 
(PDS) of black hole binaries in hard states is dominated by a component, which has a spe- 
cific shape roughly described by a broken power- law. The low-frequency part is mostly flat, 
while the power-law index a above the "break" frequency is variable between 1 and 2. 
It is well established that the fractional root-mean-square (rms) variability in a source light 
curve decreases as a source evolves from LH state to HS state. Simultaneously, both and 
a increase. Altho ugh empirical shot -noise models were able to describe in general the ob- 
served PDS shape (IFocke et al.ll2005l ). the physical picture explaining the observed evolution 
during spectral transitions was mi ssing. Moreover, shot-noise models were challenged by the 
linear absolute rms- flux relation (jUttley. McHardy fc Vaughan I l2005l : lUttley I l2004j ). This 
rms flux relation assumes that amplitudes and time-scales of shots are not independent, but 
are related in some way. 



Lyubarskii I (119971 ). hereafter L97, was the first to suggest a model for this time variabil- 



ity production in the accretion powered X-ray sources. He considered small amplitude local 
fluctuations in the accretion rate at each radius, caused by small amplitude variations in the 
viscosity, and then studied the effect of these fluctuations on the accretion rate at the inner 
disc edge. His linear calculations show that if the characteristic time-scale of the viscosity 
variations is everywhere comparable to the viscous (inflow) time-scale, and if the amplitude 
of the variations is independent of radius, then the power spectrum of luminosity fluctuations 
is a power- law l/u. If the amphtude of the variations increases with radius, the slope of the 
power spectrum of the luminosity variations is steeper than 1. Lyubarskii pointed out that 



he ha d no physical model for the cause of such fluctuations. lUttley. McHardy &: Vaughan 



(120051 ) pointed out that rms-flux relation is naturally explained in the framework developed 
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by Lyubarskii. 

TSA07 formulated and solved the problem of local driving perturbation diffusion in a 
"disk-like" configuration (which can be either a geometric thin Keplerian accretion disk or 
Compton cloud). The problem of the diffusive propagation of the space distributed high- 
frequency perturbations is formulated as a problem in terms of the diffusion equation for the 
surface density perturbations. This equation is combined with the appropriate boundary 
conditions. The formulation of this problem and its solution are general and classical. The 
parameters of the resulting PDS, diffusion time scale of the diffusion propagation of the 
local perturbations to and the power-law index of the viscosity distribution over radius, are 
essential parameters of diffusion in a given bounded configuration. In TSA07 we call our PDS 
model for the Green's function of the bounded configuration as a white-red noise (WRN) 
and we adopt this name throughout this paper. 

The problem formulation was similar to the Lyubarskii's scheme. However, the method 
of solution is different. Lyubarskii treated the factorization of the driving term (i.e. separat- 
ing it into two parts each depending on time and radius only) by linearizing the system. In 
TSA07 the analytical solution is obtained for the case of the factorized driving sources. Then, 
using the mean value theorem we showed that the general solution is simply a convolution 
of the configuration response signal (the Green's function) and the mean driving signal in 
the configuration. Thus the resulting power spectrum of the X-ray signal, as a convolution, 
is a product of the power spectrum related to the configuration Green's function and that 
related to the perturbation sources (sources of driving oscillations) . 

The PDS of the Green's function is a white-red noise power spectrum (WRN). Specifi- 
cally, the low frequency (LP) asymptotic form of the WRN PDS, when the frequency is less 
than the inverse of diffusion timescale in the "disk-like" configuration t^-^, is characterized 
by a fiat shoulder (white noise). In other words, the LP white noise shoulder is insensitive 
to the source and viscosity distributions over radius. The high frequency (HP) asymptotic 
form of WRN is a power law with index a, which is determined by the viscosity and 
perturbation source distributions over the accretion configuration. The index a — 3/2 when 
the viscosity linearly increases with radius and the perturbation sources distribution is quasi- 
uniform. The basis of the presented power spectrum formation scenario is that the timing 
signal of the WRN PDS shape is a result of diffusive propagation of driving perturbations 
in the bounded configuration (disk or Compton cloud) in the same way as X-ray photon 
spectrum is a result of the photon diffusion (namely, upscattering of seed photons) in the 
same bounded configuration. 

The driving oscillation amplitude is assumed to be a smooth function of the radius. 
TSA07 suggested that driving fiuctuations in the configuration can be introduced by g- 
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mode driving oscillations at any given annulus. The local g-mode driving fluctuations, 
produced possibly by local Rayleigh- Taylor local instabilities, are high-frequency damped 
quasi-periodic oscillations (QPOs) which frequencies are related to the local Keplerian fre- 
quencies. As we mentioned above TSA07 formulated and solved a problem of the diffusive 
propagation of the space distributed high-frequency perturbations in the bounded configu- 
ration. Our diffusion model for PDS is a product of WRN PDS and the driving source PDS 
(Lorentzian). 

The WRN PDS is a power spectrum of the solutio n of the initial value (Cauchy) problem 



which is a linear superposition of exponential shots [see lWood et al. I (120011 )]. For example, if 
the driving perturbations are distributed according to the first eigen-function of the diffusion 
operator then the bounded medium works as a filter producing just one exponential shot as 
a result of the diffusive propagation of eigen-function distribution of the seed perturbations. 
In the general case the resulting signal is a linear superposition of exponential shots which 
are related to the appropriate eigen-functions. Furthermore, TSA0 7 demonstrate that the 



observed rms-flux relations [e.g. lUttley. McHardy fc Vaughan I (120051 )] is naturally explained 



by our diffusion model. In the framework of the linear diffusion theory the emergent pertur- 
bations are always linearly related to the driving source perturbations through a convolution 
of the configuration Green's function and source distribution. 

An important question is what our diffusion model predicts for relative contributions 
of the WRN PDS, the driving oscillation PDS in the resulting PDS and for a dependence 
of the integrated PDS power on the driving oscillation frequency. The next question is how 
this model dependence of the integrated power vs the driving oscillation frequency is related 
to the observed dependence of that. The answers to these questions are the points of the 
presented study. 

In §2 we refer to details of Cyg X-1 observations with RXTE. In §3 we outline the main 
features of the diffusion model and related formulas. In §4 we show how the model integrated 
power vs the driving oscillation frequency fits X-ray data from Cyg X-1. In §5 we present the 
inferred correlation of the Reynolds number with the driving oscillation frequency and the 
spectral state (photon index). Application of the paper results to the observed index-QPO 
frequency correlations is considered in §6. Discussion and final conclusions follow in §7. 



2. Observations 



For our analysis we used Cyg X-1 da ta from t he Pr oportional Counter Array (PCA) and 
All-Sky Monitor (ASM) onboard RXTE Swank I (119991 )]. The data are available through the 
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GSFC public archive 0. In this Paper we present the analysis of a representative subset of 
■R Xr.E observations of Cyg X-1. A re ader can find the details of data reduction and analysis 



m 



Shaposhnikov fc Titarchuk I (120061 ) and TSA07. We chose approximately 200 observations 
to cover the complete dynamical range of the source evolution from low/hard to high/soft 
state. For the presented analysis we refit PDSs with our new model and we used the results 
of our previous spectral analysis for photon index F. 

To fit a PDS we used a sum of our perturbation diffusion model and one or two 
Lorentzians to account for the Low Frequency Quasi-Periodic Oscillations (LFQPOs). For 
higher photon inicies, when the source is close to high/soft state, the contribution of the 
accretion disk variability component sometimes becomes significant. It is observed as an 
additional power law at the lower frequencies (see TSA07 for details). We fit this component 
with simple power law, when it is needed. 



3. The main features of the model and underlying assumptions 

We consider a scenario related to our model (see also TSA07) where the Compton Cloud 
( "disk-like" configuration) is located in the innermost part of the source and the Keplerian 
disk is extended frorn . the C ompton Cloud (CC) to the optical companion [see Fig. 1 in 



Titarchuk. &: Fiorito I (120041 ). hereafter TF04, for the model geometry). 



We remind the reader that the diffusion equation for the surface density perturbation 
is the same for any "disk-like" configuration for which the rotational frequency Q(R) as a 



function of radius R has a Keplerian- like profile namely oc i?^^/^ [see IWood et al. I ( 20011 ) 
for details]. The standard Keplerian disk and ADAF kind of flow (Compton Corona) are 
particular examples of these accretion "disk-like" ' configurations. 

The CC presumably contracts when the source goes to the softer state. TF04 infer 
that the CC cloud in low/hard state is about 40 Schzwarzchild radii [Rs) whereas that in 
high/soft state is about 4-5 Rs- It is obvious that the variability of the mass accretion rate 
in the cloud leads to the variability of the gravitational energy release there. 

The Earth observer sees this mass accretion rate (M) variation as a variability of X-ray 
flux coming from the source given that M regulates the supply of the soft seed photons up- 
scattered off CC hot electrons. According to our scenario in the equatorial plane the plasma 
is more dense and consequently colder than that in the CC outer parts. The soft photons 
are produced in these relatively cold equatorial layers of CC and then they are Comptonized 
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off fiot electrons of tlie CC outer part forming tlie resulting X-ray spectrum. 

Tlius one can conclude that the variability of X-ray signal at energies higher than 3 keV 
(RXTE energy range) is mostly determined by the fluctuations of the luminosity ^L^it) 
originated in the CC. Here we study the luminosity fluctuations AL^.(t) in the CC configu- 
ration. 

We assume that the mass accretion rate variations AM(0, t) is converted with efficiency 
Eeff into the variations of the X-ray luminosity, i.e. ALx{t) = eeffAM{0,t). TSA07 show 
that the fluctuations of the resulting X-ray oscillation signal ALx{t) due to the diffusion of 
the driving perturbations is 

AL,{t)= [ ip{t')Y{t-t')dt', (1) 
Jo 

i.e. a convolution of the Green's function of the bounded conflguration Y{t) (WRN) and the 
source variability function (p{t). The resulting power spectrum is 

||F.H|p = ||F,Hin|FyM|p (2) 

where F^^uj), F^i^uj), Fy{uj) are Fourier transforms of AL^(t), fit), Y(t) respectively. The 
WRN PDS ||Fy(tt;)|p is described by Eq. (64) in TSA07 (also see the asymptotic forms of 
that in Eqs. 65-66 there). 

Using the total power of the driving oscillations Pdr one can present the driving oscilla- 
tion PDS as (see Eqs. 22 and B5 in TSA07) 

\\F,ml = (3) 

where = uj^r/ (27r) is the driving oscillation frequency, f is a full width of half maximum 
(FWHM) of the Lorentzian and a constant a varies in the range between 1 and 2 depending 
on the ratio of 2vdrl^dr'- 

2 arctan(2z/dr/f dr) , , . 

a = IH . (4) 

TT 

For example a = 1 and a = 2 when 2i'iij./i'dr ^ 1 and 2udr/'^dr ^ 1 respectively. 

It is worth noting that Eqs ([HH]) have been derived in TSA07 with an assumption that 
the local driving perturbations ^{t,^) are damped quasiperiodic oscillations, namely 

$(t, = A^f{i) e^Y>{-\^drt + lUJdrt) (5) 
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where the amplitude of the driving oscillation is function of R only = R^/"^). Note the 
damping coefficients F in Eq. ([3]) and F in Eq. ([5]) are related as F = 2'kV . 

Then TSA07 apply the mean- value theorem for the integral W{R, t) of Eq. B2 in TSA07. 
Namely the mean- value theorem in its general form for the product of two continuos functions 
g{t) and f{t) states that if a function g{t) is positive, i.e. g{t) > then 

\{t)f{t)dt = f{x) f g{t)dt 



where a < x <h. Given that the product of the Green's function G{R, ^, t) and the driving 
oscillation amplitude /(^) in Eq. B2 in TSA07 is positive we obtain that 



W{R,t)=A^ / dt'exp[-T,rt'/2 + zuJar{^,)t'] / G{R,^,t - t')f{Od^ 

Jo J Rin 

where i?^ < ^* < Ro- Then TSA07 proceed with the Fourier transformation and PDS 
determination of W{R,t) and ultimately they derive equations ([Hll]). 

In Fig. [T]we show a typical example of the model fit to the data using Eq. ([2]). The 
parameters of PDS continuum for white- red noise component (WRN) ||Fy(ci;)|p (see for 
details TSA07) are the diffusion time scale t^, index of the power-law viscosity distribution 
ifj and for the driving oscillation component ||F(p(ti;)|p they are Vdr and Vdr (see Eq. |3]). 

The power-law index of the viscosity distribution ip is related to the power-law index of 
the red noise in the WRN PDS (see TSA07): 

a = = —r, TT for w > U 

and 

a = 2 for ^ < 0. (6) 



The model predicted integrated power is (see TSA07) 

7-. PdT 



x,diff 



2TiPdMQ + l/AQ)DC 



(7) 



where Q = tdr/ ^dr is a quality factor for the driving signal and D is a factor of order of unity. 
Equation (j7]) was derived use the mean value theorem for the integral of the product of two 
functions. In TSA07 we assumed that a constant C related to the mean value of ||-F<^(i^)||^ 
over the frequency integration range is about a few (see Appendix B2 in TSA07 for details). 
Here we specify and obtain C— constant when we compare the model dependence Px,diff on 
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Udr with the observable Px on Udr- In order to make this comparison one should determine 
the best-fit parameter to and Pdr as functions of Udr- Note that to derive Eq. ([7]) we also use 
the fact that the WRN PDS \ \Fy{uj)\\'^ is normalized to l/iDto) where D>1. (see Eq. B16 
in TSA07). 

Now we follow the method suggested in TSA07 to infer Pdr vs Udr from the observations. 
Namely, given the fact that the driving PDS is a constant at frequencies u <^ Vdr we have 



d + (r../2)2 



Because for any power spectrum ||F(ct;)|p 

||F(cu)||2rfcu = ||F(27rz/)||'27rt^z/ = \ \F{u)\\ldu 
we obtain that (compare with Eq. [2]) 

||F.(z.)||^ = (27r)-i||F,(z.)||^||Fy(z.)||^. (9) 

Thus a combination of Eqs. (IHl) and ([9]) leads us to determination of the integrated 
power of the driving oscillations 



^<ir = ||F,(0)||^ = 27r^ali^^-41i^^H^- (10) 



+ (rd./2)2 






\Mo)W 



anjuj + {T dr/2Y^ 

We remind a reader that the values of Udr and tdr are the best-fit PDS parameters, ||Fi,.(0)|p 
is the observed PDS value at z/ = and ||Fy(0))|p is a value of the normalized WRN PDS 
at z/ = 0. As we mention above the integral of the normalized WRN PDS ||Fy(z/))|p over u 
is l/{Dto). 



4. The integrated power vs driving oscillation frequency and photon index 

In Figure [2] we present the observed correlation of low frequency QPO centroid ul with 
the driving QPO frequency Udr (upper panel) and photon index F with Udr (lower panel). 
These correlations imply that z/^ along with Udr increase when the source becomes softer. In 
other words, the emission area [Compton cloud (CC)] contracts when the source evolves to 
the soft states. 

In Figure [3] we also see this effect of CC contraction as anticorrelation of the CC diffusion 
time scale to with Udr- The inferred dependence of to vs Udr can be fitted by the power law 
to OC z/^^, 
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We infer the integrated power of the driving oscillations Pdr vs Vdr (see Eq. [TO]) and 
then we obtain the model integrated power Px,diff vs Udr (see Eq. [7]). In Figure H] we show 
that the dependence of Pdr vs Udr can be fitted by power law Pdr oc J^dr'^'^'^'^^- Namely 
the driving oscillation power Pdr decreases when the source (Cyg X-1) goes to softer states. 
Presumably the decay of Pdr with Udr is also related to the contraction of Compton cloud. 
The driving osci l lation s can result frorn the R ayleigh- Taylor (RT) local instability [see e.g. 
Chandrasekhar I (jl96ll ) and iTitarchuk I (120031 )] ]. The decay of Pdr can be considered as a 
cumulative effect of the local Rayleigh- Taylor (RT) instability when the effective area of a 
given configuration (CC) undergoing RT oscillations contracts. 

In Figure [5] we present a comparison of the observable PDS integrated power P^ (black 
filled circle) with the model predicted Px,diff (crosses) (see Eq. [7]). One can see that the 
dependence Px,diff on Udr is similar to the observable correlation P^ vs Udr- 

Note that we obtain the factor C ~ 4 (see Eq. [7j) by shifting a set of the values of Px,diff 
along Y-axis to fall on top of Px values. The power-law Px oc fits the dependence 

of the theoretical and observable integrated powers vs the driving oscillation frequency. 



5. The Reynolds number of the accretion flow in Compton cloud configuration 



Titarchuk. Lapidus fc Muslimov I (119981 ). hereafter TLM98 introduced the Reynolds num- 
ber Re = VrR/z/ (7 in their notation) where Vr, D are an average radial velocity, an 
average viscosity over a given configuration respectively and i? is a configuration scale. 
They demonstrate that the size of the transition layer (CC) between the fast rotating ac- 
cretion disk and the relatively slow rotating central object (either BH or NS) strongly de- 
pends on the Reynolds number Re. It is worth noting that a— viscosity parameter intro- 



duced by I Shakura fc Sunyaev I (Il973l ) is related to the TL (CC) parameter Re. Given that 
a ~ O/eVgH), where Vg is the sound velocity in the "disk- like" configuration and if is a half 
of the configuration vertical size, we obtain that 

Re ~ A/a. (11) 

To infer this formula one should assume that A = H/R and Vs ~ Vr. 

TSA07 has already shown that there is a possibility to infer Re, and consequently 
a— parameter (see Eq. [TT]), using the RXTE observations of BHs and NSs. This deter- 
mination can be done if a particular RXTE power spectrum (PDS) has the white-red noise 
(WRN) component and also if the QPO low frequency ul is present there. TSA07 argue that 
the relation between observed frequencies ul and 1/to seen in Cyg X-1 PDSs is similar to the 
theoretical relation between the diffusion frequency i^diff and the frequency of CC volume 
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(magneto-acoustic) oscillations vma [see e.g. iTitarchuk. fc Osherovich I (Il999l )]. This fact 
leads them to conclude that z/^ can be con sidered as the frequency of the magneto-acoustic 
oscillations which theory was developed by ITitarchuk et al. I (120011 ). hereafter TBWOl. 



Using the best-fit parameters of the PDSs we can even infer the evolution of the Reynolds 
number of the accretion flow in the Compton cloud (CC), with the change of photon index 
r. In fact, TSA07 relate diffusion time to with Re and a magneto-acoustic QPO frequency 

I'M A as 



4 



3 (4 - 



VrRq 



Rq 



4 4 



Re 



3 (4 - ?/')2 auAl^MA 



where Vn ~ K and Oma is a numerical coefficient. Formula ( |T2l) leads to the equation 

Re = 2. f ^) ^^i^(..to) 
V 27r / 4 4 ^ ^ 



(12) 



(13) 



that allows us to infer a value of Re using the best-fit model parameters to and the QPO low 
frequency pl which is presumably close to vma- Ultimately we can find the evolution of Re 
with the photon index F given that vl, to and the viscosity index %p evolve with F (see Figs. 
HI and E]). 

In Figure [7] we present the inferred Reynolds number as a function of the photon index 
F. We use Eq. ( 1T3|) where we set uma = 27r (see details of this determination of a^M in 
TBWOl and TSA07) and the observable correlations of and to with F (see Figs. [2] and 
E]). One can see that Re— number steadily increases from 10 to 70 when the source evolves 
from low/hard state to high/soft state. In contrast, TSA07 found that Re ~ 8 ± 2.5. We 
note, however, that in TSA07 only WRN model was used without accounting for driving 
oscillation distribution, which significantly affects the resulting value for ip. They also used 
a limited set of data. 

Note the observed behavior of the Re-number vs F and mass accretion rate was pre- 
dicted by TLM98, where they formulated a transition layer model (TLM) and studied its 
consequences for observations. It is important to emphasize that the Re— number along 
with the photon index F, the low frequency QPO and the driving frequency v^r can be 
considered as characteristics of the spectral state. All of them correlate with each other. 



6. Photon index-QPO frequency correlation 

TLM98 showed that the outer (adjustment) radius of the transition layer (CC) Rout 
measured in the dimensionless units with respect to Schwarzschild radius Rs = 2GM/c^, 
Tout = Rout/ Rs^ anticorrelates with Re— number or photon index F (spectral state) only. Thus 
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vl (or vma) as a ratio ~ Vma/ Rout should correlate with F (or Re— number) where values vl 
related to the same F for different sources should be inversely proportional to a mass of the 
central object (black hole or neutron star). Note that a plasma velocity Vma is also a function 
of F only. The comparison of the observed index-QPO frequency correlations for two differ- 
ent sources (with two different masses) should lead to determination of their relative masses 
with respect each other. This is the main idea behind the method of weighing black holes 
(TLM98) recently applied for BH mass determi n ation i n a number of Galactic and extra- 



galact i c sources [see TF04 ; iFiorito fc Titarchuk I (l2004f): I Dewangaii. Titarchuk fc Griffiths 



feooeh: Istrohmaver et al.~l J2007h and Ishaposhnikov fc Titarchuk I j2007lVI 



Discussion and Conclusions 



Chakrabarti fc Titarchuk I (119951 ) proposed a model of the truncated Keplerian disk and 



hot corona locate d in the innermos t part of the source using arguments in terms of radiation 



hydrodynamics. IChurazov et al. I (119991 ). hereafter CGR99, were the first who supported 
this model using the specific form of the power spectra of X-ray radiation. Particularly they 
noted that two-component X-ray spectra (soft multicolour black-body and harder power 
law) are frequently observed from accreting black holes. These components are presumably 
associated with the different parts of the accretion flow (optically thick and optically thin 
respectively) in the vicinity of the compact source. They also emphasized that for Cygnus 
X-1, the overall shape of the power density spectra in the soft and hard spectral states can 
be qualitatively explained if the geometrically thin disc is sandwiched by the geometrically 
thick corona extending in a radial direction up to a large distance from the compact object. 
In the hard state the thin disc is truncated at some distance from the black hole followed by 
the geometrically thick flow. The break in the PDS is then associated with the characteristic 
frequencies in the accretion flow at the thin disc truncation radius. 

However there is a difference between CGR99's suggestions and the results of our paper. 
We come to the model of the truncated disk plus corona using the fit of the theoretical model 
of the power spectrum (see TSA07) to the data. Whereas CGR99 proposed a qualitative 
explanation of the data using Lyubarskii's idea on the formation of the power spectrum in 
X-ray binaries [see L97]. 

Our results confirm the CGR99 conclusion that a stable optically thick (but geomet- 
rically thin) disk extends down to small radii (as indicated by the strong and stable soft 
component) and whereas prominent variations of the harder component (presumably related 
to Compton cloud) are present in the broad range of time-scales (up to at least 10^ s; see Figs. 
2-3 in TSA07). We also confirm the observed emergent spectrum consisting of two compo- 
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nents (soft stable component due to the disc emission and harder highly variable component 
due to Comptonization in the corona). The relative contribution of these two components 
to the luminosity of an averaged spectrum would then reflect the ratio of the energy releases 
in the disc and corona (or mass accretion rates). It is also worth noting that our presented 
analysis indicates that the high-frequency part of the PDS (QPO and turnover) at about 10 
Hz and higher is possibly related to the local instabilit ies operating in the Comptoniz ation 



region of the energy release (see CGR99). Note that iGierlinski fc Zdzdiarski I (120051 ) also 
pointed out the evolution of PDSs in XTE Jl 5550-564 and XTE J1650-500. In particu- 
lar they showed a variable power in the Comptonized component (the PDS high-frequency 
component) when a given source evolves from hard to soft states. 

However there is a difference in details between the CGR99 and results of our papers 
(see also TSA07). For example the break in the PDS frequency range between 0.01 and 1 
Hz is associated with the characteristic diffusion frequencies of the Compton cloud but not 
as that in CGR99 "with the characteristic frequencies in the accretion flow at the thin disc 
truncation radius." 



Axelsson et al. I (120051 ). hereafter ABL05, suggested a model of the PDS of Cyg X- 
1 which successfully presented a natural transition from hard state through intermediate 
state to soft state and back, and also allowed them to study the behavior of the Lorentzian 
components in detail. Using this model consisting of a power-law and two Lorentzian profiles, 
ABL05 were able to fit the Cyg X-1 PDSs, and to follow the components from hard state 
through the transitions and back. By choosing this simple empirical approach for the PDS 
data analysis ABS05 could model the major features in all states with only three components. 
The parameters used in ABL05 were the width and peak frequency of the Lorentzians, along 
with the power at the peak frequency. The parameter evolution and its correlation between 
each other can all be described by continuous functions. 

However ABL05 did not address to the nature of the components and continuum of the 
power spectrum. Their PDS description was pure phenomenological (compare with that in 
TSA07 and this paper). 

ABL05 also showed that the PDS of Cyg X-1 is dominated by the same two Lorentzian 
components at all times. They came to this conclusion based on the fact that the relation of 
the peak frequencies of these Lorentzian components follows the pattern seen in both other 
BH and NS systems, and therefore they proposed that the physical processes responsible for 
them cannot be explained by invoking magnetic fields, a solid surface, or an event horizon. 

It is worth noting that we also find that the diffusion (break) and QPO frequency follows 
the BH-NS pattern (see details in ST06 and TSA07). In addition to this we demonstrate that 



-13- 



there is a certain correlation between QPO low frequency and driving (high) frequency and 
photon index (see Fig. 2) and show a correlation between the diffusion time scale (an inverse 
of the diffusion frequency) and the driving frequency (see Fig. 3). Furthermore we explain 
the nature of the observational appearances of the power and photon spectra evolution by 
the evolution of the truncated disk and corona configurations when a X-ray source evolves 
from hard state to soft states. We come to the conclusion that the truncated disk and corona 
scenario should be a common model in BH and NS systems. 

In fact, TLM98 showed that the bounded configuration (Compton cloud) surrounding 
compact objects is the transition layer (TL) that is formed as a result of dynamical adjust- 
ments of a Keplerian disk to the innermost sub- Kcplcrian boundary conditions. They argued 
that this type of adjustment is a generic feature of the Keplerian flow in the presence of the 
sub-Keplerian boundary conditions near the central object and that it does not necessarily 
require the presence or absence of a hard surface. TLM98 concluded that an isothermal 
sub-Keplerian TL l)otween the NS surface and its last Keplerian orbit forms as a result of 
this adjustment. The TL model is general and is applicable to both NS and BH systems. 

As the conclusions we want to single out 

1. That we explain the decay of the emergent time variability of X-ray emission in 
compact sources when these sources evolve from low/hard to high/soft states. We find that 
the resulting power from Cyg X-1 decays with the driving oscillation frequency v^ir as 

2. We show that the reciprocal of the diffusion time scale of the perturbation ^, the 
low frequency QPO vl-i the driving oscillation frequency Vdr-i inferred by fitting the Cyg X-1 
PDSs with our diffusion model, increase when the source evolves from low/hard state to 
high/soft state. This behavior of the PDS characteristics impfies that the Compton cloud 
contracts towards softer spectral states. The driving oscillations are probably caused by the 
local Rayleigh- Taylor instability which cumulative Pdr decreases when the effective area of 
the configuration producing the RT oscillations contracts. The decay in driving power leads 
to the decay in the total observed variability power from the source. Using the fact that ^, 
^L-i ^dr increase with F and Pdri Px decrease with F we conclude, as a result of our analysis, 
that the Compton Corona shrinks when Cyg X-1 goes from low/hard state to high/soft state. 

3. Our extensive data analysis of the power spectra from Cyg X-1 with an application 
of the method of Re— number determination developed in TSA07 indicates that Re related 
to Compton cloud configuration increases from values about 10 in low/hard state to that 
about 70 in high/soft state. We confirm the predictions by TLM98 that Re— number should 
increase with index and QPO frequencies. Thus one can conclude that the observable index- 
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QPO correlation is probably driven by the increase of Re— number when the source evolves 
from low/hard state to high/soft state. It is worth noting that inverse proportionality of the 
low-frequency QPO with respect BH mass in the index-QPO correlation leads to the method 
weighing BHs employing this index-QPO correlation. 

We acknowledge very useful comments and suggestions by the referee. 
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Fig. 1. — A particular example of observable PDS. The PDS continuum is fitted by our 
diffusion PDS model which is a product of WRN PDS and the driving oscillation Lorentzian. 
We also use a simple Lorentzian to fit QPO features. 
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Fig. 3. — The best-fit model parameter, diffusion time scale to vs u^r- The dashed fine is the 
best fit power law to oc i^ar ^^^^^"^- 



-19- 




J I I I I L 



20 



50 



100 



200 



Vdp Hz 



Fig. 4. — Inferred P^r vs Udr- The dashed hne is the best fit power law Pdr oc 
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Fig. 5. — Comparison of the model Px,diff vs Udr (crosses) (using Eq. [7j) with the observable 
Px vs Udr (black filled circle) The dashed line is the best-fit power law Px oc z/^°'^^^° °^. 
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